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We investigate multicritical phenomena in 0(A'’)+0(M)-models by means of nonpertnrbative 
renormalization group equations. This constitutes an elementary building block for the study of 
competing orders in a variety of physical systems. To identify possible multicritical points in phase 
diagrams with two ordered phases, we compute the stability of isotropic and decoupled fixed point 
solutions from scaling potentials of single-held models. We verify the validity of Aharony’s scaling 
relation within the scale-dependent derivative expansion of the effective average action. We discuss 
implications for the analysis of multicritical phenomena with truncated flow equations. These hnd- 
ings are an important step towards studies of competing orders and multicritical quantum phase 
transitions within the framework of Functional Renormalization. 
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I. INTRODUCTION 

The antinomy of two equally justified principles makes 
the essence of tragedy in literature (Hegel). In the same 
way, the competition of ordering principles in complex 
systems adds a new enthralling element to their macro¬ 
scopic phenomenology. The purest manifestation of this 
interplay is two either intertwining or mutually excluding 
ordered phases in a physical system mu- Correspond¬ 
ing examples can be found everywhere in physics, includ¬ 
ing magnetic systems [TMT^ . high-temperature super¬ 
conductors [I7H22], graphene [23U26] . and dense quark 
matter [?7l - l^ . The theoretical description of such sit¬ 
uations is complicated by the need to accurately resolve 
both universal and non-universal features of the particu¬ 
lar phase structure. 

In fact, many aspects of competing orders are already 
universally dictated by the theory of critical phenom¬ 
ena. For this consider a physical system with two dis¬ 
tinct macroscopically ordered phases which are separated 
from the disordered phase by means of second order phase 
transitions each. The corresponding phase transitions are 
approached by separately fine-tuning two parameters, gi 
and 32 , which can be temperature, coupling strength, 
etc. One may now, for instance, want to know whether 
there is a second order multicritical point ( 31 c, 32 c) in 
the phase diagram where both transition lines meet, or 
whether there is a coexistence phase of both orders. 

If a second order multicritical point exists, the sys¬ 
tem in its vicinity will be described by a critical quan¬ 
tum field theory for the two order parameter fields. The 
properties of this field theory also determine whether 
there can be a coexistence phase. Whereas classical mul¬ 
ticritical phenomena are typically captured by an eu¬ 
clidean 0(A^)-|-0(M)-model, bosonic self-energy effects 
or the presence of gapless fermions can complicate the 
setting for quantum phase transitions [50115^ . Further¬ 
more, even if a multicritical point exists, the microscopic 
parameters of the model may not lie in the basin of at¬ 
traction of this fixed point, or there may be energetic 
reason for the second order lines to actually become first 


order lines before they meet. In these cases a coexistence 
phase is excluded as well. 

Thus the theoretical modelling of competing orders 
requires both a resolution of universal properties of a 
given system under consideration, e.g. critical phenom¬ 
ena, but also non-universal ones, e.g. the effective poten¬ 
tial to resolve first order phase transitions. Furthermore, 
the method should be applicable beyond the realm of 
classical phase transitions. A versatile tool which quali¬ 
fies here is provided by the Functional Renormalization 
Group (FRG) [55M5] . It naturally captures universal 
aspects and critical phenomena in the long-wavelength 
limit, but at the same time allows to resolve system- 
specific properties at intermediate scales. In the fol¬ 
lowing we focus on classical multicritical phenomena in 
0(iV)-l-0(M)-models. Those are the key building blocks 
for more involved setups as the bosonic two-field model 
always appears as a subsector of the corresponding set of 
beta functions. 

In the vicinity of one of the second order lines de¬ 
scribed above, say the one associated to tuning 31 , the 
corresponding order parameter fluctuations are captured 
by an effective action or free energy functional F for 
the order parameter field 3>, which we assume to be 
0(A^)-symmetric. We approximate here the effective ac¬ 
tion within the local potential approximation including 
a wave-function renormalization (called LPA' in the fol¬ 
lowing) and write 

r^[$] ~y’d^a:(-Z$<i>-U($)), (1) 

where $ is an iV-component vector and <& = |$|. U(‘I>) is 
called the effective potential. The translation invariant 
ansatz is also applicable in the low-energy limit for spin 
systems close to the phase transition. The critical point 
of the effective action in Eq. ([^ is approached by fine- 
tuning a parameter ct —^ CTc- Close to the transition line 
we have a linear relation gi — (x a — ac which links the 
actual physical system to the corresponding 0(A^)-model. 
In the same way, the vicinity of the multicritical point 
is captured by means of an 0(A^)-|-0(M)-model, which 
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we approximate within LPA' according to the effective 
action 

+ !/($, d/)). (2) 

Herein, $ and d' are N- and M-component vectors, 
respectively. Furthermore we restrict to the three- 
dimensional case in the following, d = 3. 

A reasonable candidate for a second order multicriti- 
cal point of the two-field model is given by the decoupled 
fixed point (DFP) solution. In this case = 

F^[<h]-l-r'^['h], such that the critical properties are inher¬ 
ited from the individual single-field models. This solution 
represents a stable fixed point of the theory provided that 
mixed terms in the effective action, like f 
become irrelevant in the infrared. In particular, the scal¬ 
ing dimension of the coupling A$^, denoted by 63 in the 
following, is negative for a stable DFP. Aharony’s exact 
scaling relation [51 states that 6*3 satisfies 

= 1 + 1 - d (3) 

1/1 V 2 

at the DFP. Herein, and 1^2 are the usual correlation 
length critical exponents of the associated 0{N)- and 
0(M)-models in d dimensions, respectively. The latter 
are well-studied in the literature | 6 ] , from what the possi¬ 
bility of a stable DFP can be deduced immediately. Note 
that the DFP solution is such that it allows for a coexis¬ 
tence phase of both orders em¬ 
it turns out that for the important cases N,M < 3 
the value of 63 is rather small. Thus, quantitatively 
small errors in the Vi or violations of the scaling re¬ 
lation in the two-field model might turn a stable fixed 
point into a seemingly unstable one. Previous studies of 
0(A^)-|-0(M)- and 0 (A^i)-|- 0 (A 2 )-l- 0 (A^ 3 )-models with 
the FRG within LPA' indicate a small but visible viola¬ 
tion of the scaling relation [iniiiiiiim- We resolve this 
issue and show that it is related to a truncation-related 
ambiguity in defining the stability matrix of the critical 
theory. We demonstrate that the scaling relation is valid 
within LPA' for a commonly used regularization scheme. 
We argue that although the scaling relation might be vi¬ 
olated, the correct fixed point is approached during the 
renormalization group flow. This is of great importance 
for computing the phase structure of more complicated 
models. 

This paper is organized as follow. In Sec. jlTjwe discuss 
critical phenomena of 0(7V)-models in three dimensions 
by determining the associated scaling potentials. Those 
are then used in Sec. Ell to study the stability of the 
decoupled and the isotropic fixed point in 0(iV)-|-0(M)- 
models. In particular, we compute the critical expo¬ 
nents of the associated multicritical points and show that 
Aharony’s scaling relation is valid within LPA'. In Sec. 
|IV| we comment on an ambiguity in defining the stability 


matrix for the truncated system, and show how it in¬ 
fluences both critical and multicritical phenomena. We 
summarize our main findings in Sec. jVjand give an out¬ 
look on possible extension of our approach. In App. 
we derive explicit algebraic expressions for the beta func¬ 
tions of the single-field model which are used throughout 
the work. Some formulas for the two-field model which 
are used in Sec. |lV]are collected in App. 

II. SCALING POTENTIALS FOR 0(N)-M0DELS 

We compute the critical scaling effective actions F^ 
and in Eqs. Q and (j^ by means of the FRG. 

The latter is formulated in terms of the effective average 
action, F^, where fc is a momentum scale. The effective 
average action interpolates between the microscopic ac¬ 
tion on large scales and F= T |3H |35] . Its evolution 
with k is given by the exact Wetterich flow equation 

9/cFfc = iTr(-7^y^*^- dkRk]- (4) 

For an introduction to the method we refer to Refs. (31- 
l43] . Here we note that the second functional derivative 
F^ ' appearing on the right hand side makes the equation 
highly coupled and non-linear. Furthermore, the setting 
requires to specify a regulator Rk, which regularizes the 
infrared properties of the functional integral. Here we 
employ the so-called optimized regulator for the individ¬ 
ual 0(A)-fields jUlllSj. It is diagonal in field space and 
given in momentum space by 

Rk{cf) = Z^{k^-<f)Q{k^-cf), ( 5 ) 

where 0 is the step function and q is the euclidean mo¬ 
mentum. 

We seek scaling solutions of the Wetterich equation, 
which describe the system at its critical point. We intro¬ 
duce our notation for the single-field model here, but all 
notions immediately generalize in a straightforward way 
to the two-field model. By defining the dimensionless 
renormalized field 6 = —the dimensionless 

effective potential 

Vk{cl>) = ( 6 ) 

is a function oi (j) = |(/)| due to 0(A^)-symmetry. The scal¬ 
ing solution satisfies v{(j)) = 0 for all (j), where the 
dot denotes a derivative with respect to renormalization 
group time t — log(fc), i.e. Vk = kdkVk- The scale depen¬ 
dence of the wave function renormalization is encoded in 
the anomalous dimension 77 . We have 

V = ( 7 ) 

In the following we resolve the functional form of the ef¬ 
fective potential u( 0 ), but approximate the (/)-dependence 
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of r] by evaluating it at the minimum of n((/)), see Eq. 

This constitutes a truncation of the more general 
function r]{(j)). With the described parametrization and 
regularization scheme we truncate the exact equation 
such that exact relations will typically be violated. On 
the other hand, this allows for an approximate solution 
of the flow equation. 

The scaling potential for the 0(A^)-model within LPA' 
from the ERG with the optimized regulator is found as 
the solution of 


0 = + acl)v\4>) + 


b{N - 1) 

( 8 ) 


where 


a = 


d — 2 + 7] 
2 


b = 


1 \ 

d V d + 2)’ 


1 

- p+Tp72fp72)’ 


(9) 

( 10 ) 


and the dimension is set by d = 3. The value of 
the anomalous dimension rj has to be determined self- 
consistently such as to satisfy 


8vd 

d k2(1 + A)^’ 


( 11 ) 


where k is the minimum of and A = v"{k). The 
boundary conditions for the potential at the origin are 
n'(0) = 0 and w"(0) = a. Herein, cr is a relevant pa¬ 
rameter which has to be fine-tuned such that there exists 
a solution of Eq. Q for all </>. Of course, i;"(0) = a 
translates to r'(O) = bN/d{l + a). One could subtract a 
suitable term from Eq. Q such that t!(0) =0 since n(0) 
has no physical significance. 

Eq. ([^ has a trivial constant solution with rj = a = 0, 
which corresponds to the Gaussian fixed point. The con¬ 
stant solution, however, is unstable with respect to per¬ 
turbations due to the operators (jP and (/>^. The number 
of relevant (dangerous) perturbations of the solution de¬ 
cides on the likeliness of the latter to be realized in a 
physical system. We seek solutions to Eq. § which 
have only one relevant direction. In three dimensions, 
such a solution exists for all TV. It corresponds to the 
Wilson-Eisher fixed point of the renormalization group 
trajectory. To simplify the following discussion, we re¬ 
fer to the “solution” of Eq. (|^ as the one solution z;*(^i)) 
which has only one relevant direction. The latter requires 
fine-tuning of a. Further solutions might exist, but are 
not of interest to us here. 

For large field amplitudes the solution of Eq. (|^ be¬ 
haves as v{(f)) ^ . This scaling regime, however, can¬ 

not persist to small values of 4 > as the term 1/(1 + v”{4>)) 
breaks the invariance of w(/>) with respect to a rescaling 
of 4>. Furthermore, the minimum k turns out to be such 
that K < 0.5 for all relevant cases. Hence it is dominated 
by the region of small field amplitudes, and, accordingly, 
also A and ij can be deduced from a sufficiently precise 
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FIG. 1: Scaling potential for the Ising model {N — 1) deter¬ 
mined from Eq. (|^. The solid (blue) line shows «*(<(>) ob¬ 
tained with the shooting method. The dashed (red) lines dis¬ 
play the corresponding Taylor expansion with CTc found from 
the shooting solution. We show expansions to order tfP, (fp, 
and 0®. The polynomial ansatz to order (fP cannot be distin¬ 
guished from the solid line within the resolution of this plot. 
The inset shows approaching a constant value for 

large field amplitudes. 


resolution of the function for small (f). This obser¬ 
vation provides the basis for Taylor expansion schemes 
based on the shooting method [491452] . 

For the shooting method, Eq. (§ is treated as the evo¬ 
lution of v{(j)) in the formal time-variable with initial 
conditions w'(0) = 0 and n"(0) = a. For a detailed intro¬ 
duction to the method we refer to Ref. m and restrict 
to a simple sketch here. Given a value for the parameter 
r], say ry = 0, there exists exactly one a = < 7 ( 77 ) < 0 such 
that Eq. Q can be integrated up to very large (/, possi¬ 
bly (j) ^ 00 . In practice, the right a(r]) can be found to 
arbitrary precision by scanning possible candidates for a 
in nested intervals. The obtained potential ?;((/), ry) will 
typically not satisfy Eq. 0 for the anomalous dimen¬ 
sion. However, by iterating this step while using ry from 
the previous step, the solution collapses to the scaling so¬ 
lution Vi,{(j)) = v{4>, rjP) rather quick in three dimensions. 

The scaling potentials for 0(A^)-models found from the 
shooting method are numerically exact within LPA'. We 
show an example for = 1 in Fig. However, the im¬ 
portant region around the minimum k, which determines 
ry, is also well-captured by a Taylor expansion 

( 12 ) 

n=0 


around (j) = 0 with sufficiently large L. The coefficients 
Vn can be obtained from inserting this ansatz into Eq. 
Q with the right rj and (Tc known from shooting. In 
fact, the coefficient can be expressed as an explicit 
polynomial in the coefficients u„_i,..., uq, see Eq. (A7|. 
As a consequence, the full set of coefficients {vn} can 
recursively be determined from V 2 = Uc, which allows for 
very large im- Our procedure for a given 0(Al)-model 
is thus as follows: 
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( 1 ) 

( 2 ) 


Determine Uc and ry with the shooting method. 


Compute the coefficients from the algebraic re¬ 
cursion relation (A71 starting from = 0 and 

V2 = CTc- 


Model 

Cc 

V 

01 


0 ( 1 ) 

-0.16574 

0.0443 

1.545 

0.647 

0 ( 2 ) 

-0.20460 

0.0437 

1.435 

0.697 

0(3) 

-0.23588 

0.0409 

1.347 

0.742 


The explicit recursion formula (A7) for the coefficients 
Vn allows for a fully algebraic solution of Eq. ^ as 
well [S^. We name this Taylor shooting method. In 
this approach we determine a (if) such that the (squared) 
right hand side of Eq. Q, after inserting the polynomial 
ansatz ( 12 ) with a given ct, is smaller than a certain £ > 0 
for all (j) < fmax- For large enough L, the result for a 
converges to Cc as obtained from the shooting method. 
Accordingly, both methods yield the same Taylor expan¬ 
sion (12 1 and the same critical exponents. Although the 


Taylor shooting method is conceptually interesting, as 
it is fully algebraic, the shooting method is much more 
efficient in finding for practical purposes. 

The flow equation for the coefficient Vn is given by 
Pn = ( 9 ”/ 3 / 9 (/>") 0 =O) where /3((/)) is the right hand side 
of Eq. ([^ . The same algebraic recursion relations which 
dete rmine the set {u„} can be used to find {fn}, see Eq. 
( |A5[ ). We define the stability matrix Ai of the set of 
differential equations for {vn} by its entries 




dP„ 


dv„ 


(13) 


The derivative is applied for fixed rj. Due to the overall 
minus sign, relevant infrared directions are signalled by 
positive eigenvalues of M. We order the set of eigen¬ 
values {9i} of Ai such that 6 i > 82 >■■■> 9l. For 
sufficiently large L we have 


V = 


1 

Ti’ 


(14) 


where v is the usual correlation length exponent. We dis¬ 
play the values of 77 = 77 * and v found from the shooting 
method in Tab. |Tj For comparison we display reference 
values from Monte Carlo simulations and the £-expansion 
with £ = 4 — d in Tab. |llj The Gaussian fixed point with 
v((j)) = vq yields = 1/2 and 77 = 0. In this case Ai has 
eigenvalues ( 2 , 1 , 0 , — 1 , — 2 ,...). 


III. STABILITY OF FIXED POINTS FOR 
0(N)+0(M)-M0DELS 

The scaling solutions found for 0(N)-models can be 
used to obtain information on multicritical phenomena in 
0(N)-|-0(M)-models. More precisely, we can deduce the 
stability of the isotropic fixed point (IFP) and decoupled 
fixed point (DFP) for the two-field model. The input 
for this consists in the values of ac(N) and rj(N) from 
the 0(Af)-models. They allow to construct the of 
the 0(A^)-model scaling solution u*((/)), and from this we 
compute the stability matrix of the 0 (Y)-|- 0 (M)-model 
at the IFP and DFP, respectively. 


TABLE I: Critical exponents for three-dimensional 0(A)- 
models within LPA' from the FRG with the optimized reg¬ 
ulator. The stability matrix eigenvalue 9i and the critical 
exponent z/ are evaluated from the stability matrix of a poly¬ 
nomial ansatz to order L = 30 in Eq. ( |12[ |. The results are 
well-converged at this order. The value of CTc can be computed 
efficiently with the shooting method. For demonstration pur¬ 
poses we show the leading digits of CTc, although we take into 
account more digits for the numerical evaluation. 



Monte Carlo 

e-Expansion 

Model 

P 

V 

V 

f] 

0 ( 1 ) 

0 ( 2 ) 

0(3) 

0.63002(10) 

0.6717(1) 

0.7112(5) 

0.03627(10) 

0.0381(2) 

0.0375(5) 

0.6290(25) 

0.6680(35) 

0.7045(55) 

0.0360(50) 

0.0380(50) 

0.0375(45) 


TABLE II: Reference values for the critical exponents ly and 
77 of the 0(A)-model from Monte Carlo simulations and the 
£-expansion with e = 4 — d. The Monte Carlo values for 
A = 1, 2, 3 are taken from Refs. [M| . [5^ . [5S]. The results of 
the e-expansion are from Ref. [53 (labelled as “free” therein). 
The deviation of the reference critical exponents from the ones 
obtained in this work are not related to the shooting method, 
but root in the simplified LPA'-ansatz for the scaling effective 
action in Eq. 0- 


In the following we consider the two-field model along 
the lines of the single-field model, with an analogous no¬ 
tation. For the two dimensionless renormalized fields we 
write if and x to distinguish them from the field (f of the 
single field model. We also label them by “1” and “2”. 

The beta function for the effective potential v{if,x) of 
the 0(A)-|-0(M)-model from the FRG within LPA' [10] 
is given by 


P{'f,x)= - dv + aiijjv^ + a2XVx 

_l_ ^1(1 + 'I’xx) + ^2(1 + Wj/Ji/)) 

(l + 7 ;^V')(l + %x) - '''lx 

6i(A-l) b2iM-l) 


(15) 


where a subscript denotes a partial derivative with re¬ 
spect to the corresponding variable. The expressions for 
01^2 and bi 2 coincide with those for a and b in Eq. 0 
when replacing 77 —>■ 771^2 ■ 

Finding the scaling solution Vi,(if,x) from f{ip,x) = 
0 for all 7 /’,x is a rather involved task. However, two 
candidates which immediately come into mind are the 
ones associated with the IFP and the DFP. The scaling 
potential v^,(ip,x) of the two-field model at the IFP is 
given by 

^'*(V’,X) = (16) 
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where (jP = and {(/)) is the scaling potential 

of the 0{N + M)-model. The system then possesses an 
enhanced symmetry. At the DFP, instead, we have 

= + (17) 

where and are the scaling solutions for the single- 
field 0{N)- and 0(M)-models. Both field variables are 
then independent from each other. At the IFF we have 
1 i = V 2 = , whereas rji = and 772 = at the 

DFP. 

In order to decide on the stability of a scaling potential 
Vi,{ip,x) we compute the corresponding stability matrix 
M. of the two-field model. A stable fixed point has only 
two relevant directions and thus there are exactly two 
positive eigenvalues of M.. The sign of the third eigen¬ 
value, 03 , decides on the stability of x)- To compute 
the stability matrix, we insert the formal ansatz 


L L-r, 




n—0 m—0 


n\m\ 


rx^ 


(18) 


into Eq. 
cient Vnm 


(151 to obtain the beta function for the coefh- 


as 


dn Qrr. 

dij}^ dx' 


r/3(V',x) 


b=x=o 


(19) 


The entries of the stability matrix associated to the so¬ 
lution Vi,{ip,x) are given by 


M 


nm,n'm 


'b*] 






( 20 ) 


As in the single-field case, we keep 771 and 772 fixed when 
applying the derivative. 

We write f*((/i') = J2n scaling solution of 

the single-field 0(iV)-model. All Vn with odd n vanish. 
The scaling potential at the IFP with N = M then has 
coefficients 


V2n,2m 


(2n)! { 2 m)\{n-\-m)\ 2 n 
{ 2 n + 2 my.n\ml ^ 2 n-h 2 m- 


( 21 ) 


At the DFP we have 


N c- 

= v„ 5, 


mO 




nO 


( 22 ) 


with Kronecker delta 5nm- 

We summarize the results for the eigenvalues of the 
stability matrix at the IFP and DFP in Tabs. in and 
m We confirm the expectation that 62 at the IFP and 
01 and 02 at the DFP can be deduced from the knowledge 
of the eigenvalues of the critical single-field models. (The 
same is true for further irrelevant eigenvalues, which we 
do not display here, see for instance m-) The remaining 
eigenvalues, however, are genuinely determined from the 
two-field model, and are thus a consequence of the ap¬ 
proximations which lead to Eq. (151 for the beta function 
of the 0(A^)-|-0(M)-model. 


Model 

01 

02 

03 

0 (1)+0(1) 

1.756 

1.435 

-0.042 

0(1.15)-b0(1.15) 

1.767 

I. 4 O 6 

-0.001 

0(1.25)-b0(1.25) 

1.774 

1.388 

0.025 

0(1.5)-b0(1.5) 

1.790 

1.347 

0.086 


TABLE III: Stability of the IFP for the three-dimensional 
0(A)-|-0(M)-model. As the associated scaling potential only 
depends on the solution of the single field 0{N + M)-model, 
we can set N = M for simplicity. We display the three largest 
eigenvalues of the stability matrix of the two-field model at 
the IFP. The second one (in italics) coincides with 0i from 
the 0{N + M)-model. Within our approximation the IFP 
gets unstable at a critical field value of {N + M)c — 2.3, 
where O 3 becomes positive. 


Model 

01 

02 

03 

g(scal.) 

0 (1)+0(1) 

1.545 

1.545 

0.090 

0.090 

0 (l)+0(2) 

1.545 

1.435 

-0.020 

-0.020 

0(l)+0(3) 

1.545 

1.347 

-0.108 

-0.108 


TABLE IV: Stability of the DFP for the three-dimensional 
0(A)-|-0(M)-model. The eigenvalues 0; are derived from the 
stability matrix of the two-field model at the DFP. The first 
two of them (in italics) coincide with the largest eigenvalue 
of the corresponding 0{N)- and 0(M)-model, see Tab. It 
is a nontrivial finding that 03 coincides with ! = 0^ -\- 
02 — 3 from the scaling relation ©• This shows the consistent 
treatment of critical and multicritical phenomena with the 
FRG within LPA'. The scaling relation was found to be valid 
for all cases which have been tested, also including non-integer 
N and M. 


In the case of the IFP the enhanced symmetry cor¬ 
responding to the 0{N + M)-group allows to choose 
N = M without loss of generality. We find the criti¬ 
cal number of field components for the stability of the 
IFP to be given by {N + M)^ « 2.3 in agreement with 
the FRG-study in Ref. [TU]. Using an exponential regu¬ 
lator function, the upper boundary has been found to be 
(N + M)c = 3.1 with the FRG [58]. This may be com¬ 
pared to {N -b M)c = 2.89(4) and {N + M)^ = 2.87(5) 
from resummed six-loop perturbation theory and con¬ 
strained five-loop e-expansion jSHj . These results strongly 
suggest that the 0(l)-l-0(l)-model has a stable IFP. If 
N + M slightly exceeds the critical value, the stable fixed 
point is given by the biconical fixed point (BFP) instead 
(8] [To] . The BFP supports a coexistence phase and is 
linked to tetracritical behavior. The corresponding criti¬ 
cal exponents at the multicritical point differ from those 
of the IFP and DFP. 

The stability region of the DFP in the {N, M)-plane is 
bounded from below by the sign change of the third eigen¬ 
value, 03 , of the stability matrix of the two-field model. 
We find that 03 at the DFP indeed agrees with the predic- 
tion for ) from the scaling relation in Eq. M. Here 
we say that the scaling relation is satisfied if it is true 
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for the number of significant digits displayed in Tab. IV 


This accuracy is probably sufficient for deciding on the 
stability of fixed points in all relevant physical cases. In¬ 
deed, from the results presented in Tab. |^it appears to 
be very unlikely that the fixed point of an actual physical 
system is such that | 6 > 3 | < 0.001. Although the present 
analysis indicates that 6 ^ = jg valid to higher ac¬ 

curacy, an extensive numerical analysis of the associated 
convergence is outside of the scope of this work. 

Within LPA' we find that 0(V)-|-0(M)-models with 
integer values of N-\-M > 2 support a stable DFP. Due to 
the validity of the scaling relation, this information can 
also be deduced from the critical exponents Vi in Tab. 
|lj The slight deviation of the Vi from the world’s best 
values [B] due to our approximation can have profound 
implications on the (apparent) stability of the DFP in a 
given model. In particular, for the 0(l)-|-0(2)-model we 
find the DFP to be the stable fixed point. In contrast, 
from the Monte Carlo reference values in Tab. El we find 
^(scai.) ^ Q Qg 0 (l)-|- 0 ( 2 )-case, hence suggesting 

an unstable DFP. For a discussion of the current status 
of understanding of the stable fixed point in this model 
we refer to Ref. m- 

Whereas the scaling potentials of the IFP and DFP are 
constructed from the solutions of the single-field 0{N)- 
models, the determination of scaling potentials which are 
genuinely multicritical, such as the one related to the 
BFP, requires to solve the partial differential equation 
x) = 0 foi' x) in both variables, ^ and y. How¬ 
ever, inserting the Taylor expansion (18 1 into this dif¬ 
ferential equation does not yield a recursion relation for 
the Vnm- Hence the Taylor shooting method cannot be 
applied in a straightforward fashion. This also suggests 
a failure of the shooting method. The generalization of 
the shooting method to find generic scaling potentials 
x) in both field variables constitutes an interesting 
and important direction for future works on multicritical 
phenomena with the FRG. 


IV. TRUNCATED STABILITY MATRICES 


In this section we discuss a truncation-related ambi¬ 
guity in the definition of the stability matrix AI, which 
leads to an apparent violation of the scaling relation (§ 
in one case. We further show that this violation of the 
scaling relation does not conflict predictions on the sta¬ 
bility of the IFP and DFP since in both cases 6*3 turns 
out to be independent of the corresponding definition. 

For this we recall that in Eq. (131 we defined Ainm = 
—dj3nldvm by means of a derivative where the anomalous 


dimension 77 is kept fixed. However, due to Eq. 0 , 
77 depends on k and A = v"{k). With this, 77 implicitly 


depends on the coefficients Vn, and the tin-derivative may 
be applied to ri{K, A) appearing in /3n as well. In fact, 
both approaches for computing the stability matrix are 
applied in the literature. 

To make the following discussion more transparent we 


introduce u{p) = v{y/^) and the expansion 

K 

= (23) 

n—2 

where p = ((>^/2 and po = such that u'{po) = 0. The 
running couplings of the system are 

{9n} = {P0,U2,---,Uk)- (24) 


The corresponding flow equations are given by po = 
-/3'(Po)/u2, Un = /3^”Hpo) + Un+iPo- Here /3(p) is the 
right hand side of Eq. ^ expressed in terms of u. We 
have 


Pip) 


du{p) 2apu (p) + 


b 

1 -I- u'{p) + 2pu"{p) 


h[N -1) 

1 -b u'(p)' 


(25) 


The term 77„+ipo originates from the implicit dependence 
of Un on the expansion point po. The anomalous dimen¬ 
sion ( 11 ) in this parametrization reads 


167 ;^ PquI 
d ( 1 -|- 2 u 2 Po)^ 


(26) 


The coefficients Un of the scaling solution 77 * (p) are easily 
derived from the Vn of 7 ;*((^). The expressions for the 
beta functions of the 0(V)-|-0(M)-model in terms of u 
are presented in App. 

Our goal is to compare results for critical exponents 
from two different definitions of the stability matrix for 
0{N)- and 0(V)-|-0(M)-models, which we denote by 
(A) and (B) in the following. In the single-field case we 
either compute the stability matrix for fixed anomalous 
dimension, i.e. 


(A): Mnm = -(^) , (27) 

or we account for the implicit running coupling depen¬ 
dence of 77 via 


(B) : M 



(28) 


Both expressions are evaluated for 77 *. In the last line, the 
second term yields a nonvanishing contribution for pm, = 
Po or Pm = 772 . In the previous sections we have applied 
scheme (A). The generalizations of both definitions to 
the two-field model are given in Eqs. (B9) and (BIO). 

We hnd the correlation length exponents Hi obtained 
from (A) and (B) to deviate at the percent level, see 
Tab. 1^ As a consequence, we obtain different results for 




(scai.“^^^ Eq. ([^ in 0(V)-|-0(M)-models at the DFP. 
Most dramatically, applying the scaling relation yields 
contradictory results on the stability of the DFP for the 
0(I)-|-0(2)-model. 
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Model 

LPA 

LPA', (A) 

LPA' 

. (B) 


V 

V 

V 

V 

0 

V 

0 (1) 

0 

0.650 

0.0443 

0.647 

0.0443 

0.637 

0 (2) 

0 

0.708 

0.0437 

0.697 

0.0437 

0.686 

0(3) 

0 

0.761 

0.0409 

0.742 

0.0409 

0.732 

DFP 

03 

^(scal.) 

03 

g(scal.) 

03 

^(scal.) 

0 (1)+0(1) 

0.079 

0.079 

0.090 

0.090 

0.090 

0.141 

0 (l)+0(2) 

-0.048 

-0.048 

-0.020 

-0.020 

-0.020 

0.028 

0(l)+0(3) 

-0.147 

-0.147 

-0.108 

-0.108 

-0.108 

-0.063 


DFP is independent of the definition of the stability ma¬ 
trix. As a result, when considering a particular physical 
system by means of the flow equation the underly¬ 
ing 0(fV)-|-0(M)-model for low k resolves the stability 
of the IFP or DFP (i.e. the value of 63) in a unique fash¬ 
ion. The scaling relation ([^ allows to predict the value 
of 63 at the DFP by inserting the i/i of the 0(fV)-model 
as obtained from scheme (A). 

V. CONCLUSIONS AND OUTLOOK 


TABLE V: We compare critical exponents computed from 
truncated stability matrices where the variation of the anoma¬ 
lous dimension is either neglected or respec ted, l abel l ed (A ) 
and (B). They correspond to Eqs. (27l,(B9l and (28l,(B10l, 
respectively. Within scheme (A), which has been applied in 
the previous sections, and which also includes LPA as a spe¬ 
cial case, the scaling relation 63 = 9 i +62 — d is satisfied at the 
DFP. In contrast, the scaling relation is violated when apply¬ 
ing (B). However, the eigenvalue O 3 found from the two-field 
model coincides in both cases. 


In the two-field model, the ambiguity in defining A4 
seems to only afflict those eigenvalues which are inherited 
from the single-field models. Those have been highlighted 
in italics in Tabs. Ill and IV Indeed, at the IFP, the 
results for 9i and 6*3 computed with (A) agree with the 
results of Ref. m computed with (B). In contrast, the 
value of 02 , which is identical to Xjv oi the corresponding 
0(A -I- M)-model, disagrees due to the difference in the 
correlation length exponents. 

A similar behavior is found at the DFP. We find the 
first two eigenvalues 9\ and 62 within (A) and (B) to 
coincide with Xjvi from the individual single-field models 
within (A) and (B), respectively. For (A), this is shown 
in Tab. |IV[ The third exponent at the DFP, ^ 3 , is found 
to be independent of the prescription for computing the 
stability matrix. Consequently, the scaling relation 63 = 
9i + 62 — d is satisfied for (A), whereas it is violated for 
(B). Within a truncation which neglects the anomalous 
dimension, referred to as LPA, the scaling relation is also 
satisfied. These findings are summarized in Tab. 

In the analysis put forward so far we focussed on solv¬ 
ing fixed point equations. The corresponding set of equa¬ 
tions arises in the A: —> 0 limit of the renormalization 
group flow of Ffc for a system tuned to criticality. During 
the evolution with k > 0 , the apparent ambiguity be¬ 
tween (A) and (B) is lifted. The running of couplings 
does resolve the a- and A-dependence of the anoma¬ 
lous dimension. Variations in Kk and Xk will influence 
Vk = ilinkjXk) and might drive the system away from 
criticality. Hence the scaling exponents of eigenpertur- 
bations which appear during the flow are determined by 
scheme (B). Accordingly, the scaling relation is violated 
within LPA' during the flow for a Taylor expansion of 
the effective potential. 

It is an important finding of our investigation that the 
scaling dimension 03 of the operator at the IFP and 


In this work we have investigated scaling solutions 
for 0(A)- and 0(A)-|-0(M)-models in three dimensions 
within the framework of Functional Renormalization. 
The results have been obtained for a specific, but com¬ 
monly applied truncation and regularization scheme for 
the effective average action. Our main findings are sum¬ 
marized in the following list. 

(1) By a combination of shooting and algebraic recur¬ 
sion techniques it is possible to efficiently determine 
scaling solutions for 0 (A)-models and to decide on 
the stability of the IFP and DFP for multicritical 
0 (A)-|- 0 (M)-models. 

( 2 ) Aharony’s scaling relation (§ is valid at the DFP 
within LPA'. In particular, the value of 03 at the 
DFP within our truncation can be deduced from 
critical exponents of 0 (A)-models. 

(3) Previously found violations of the scaling relation 
are related to an ambiguity in defining the stability 
matrix of the truncated system. The scaling rela¬ 
tion is violated during the renormalization group 
flow of running couplings. 

(4) The value of 03 at both the IFP and DFP is not 
affected by this ambiguity. Therefore, the stability 
of multicritical points is faithfully captured by the 
running of couplings. The violation of the scaling 
relation during the flow is thus reduced to a mere 
little blemish. 

We conclude that Functional Renormalization provides 
a consistent picture of both critical and multicritical phe¬ 
nomena for scalar theories. In particular, the present 
truncation scheme with a Taylor expansion of the effec¬ 
tive potential is simple enough to be applied to more com¬ 
plicated physical systems. We do not expect the regular¬ 
ization scheme to qualitatively change any of the above 
statements. 

To reach higher quantitative precision for critical ex¬ 
ponents in 0(A)- and 0(A)-|-0(M)-models, further im¬ 
provements of the truncation are required. For one, the 
Taylor expansion of v{(j)) is expected to have a finite ra¬ 
dius of convergence. It thus fails to resolve the asymp¬ 
totic scaling behavior for large (j). Although the latter 
is captured by the shooting method, this information 
is eventually lost when defining the stability matrix in 


























terms of the associated Taylor coefficients A possi¬ 
ble way around this problem is to use spectral methods 
m, such as projection onto a complete set of orthogonal 
polynomials, to resolve the full functional form of v{(j)). 
The corresponding techniques can also be used to deter¬ 
mine the scaling dimension of eigenperturbations dv{(j)) 
of the system. It will be interesting to study whether 
approximation schemes beyond a Taylor expansion fulfill 
the scaling relation also within (B), and thus eliminate 
the above-mentioned blemish. 

Another direction of improving the present truncation 
consists in the inclusion of a field-dependent wave func¬ 
tion renormalization, or kinetial, Z$(<i>). Among other 
changes on the right hand side of Eq. (|^, we will have 
a held-dependent anomalous dimension r]((j)). The scal¬ 
ing solution then consists of two functions, and 

The corresponding field-dependence can now be 
computed with the same methods as introduced above. 
In particular, we may hope to resolve the ambiguity be¬ 
tween (A) and (B) in this way. Indeed, it probably orig¬ 
inates from the rather crude approximation « 7 ?(k) 
for all 4 >- The remaining dependence of ??(k) onto the run¬ 
ning couplings is very limited in its ways to react onto 
perturbations. Therefore, defining the stability matrix 
according to (B) need not be a consistent improvement 
of the truncation, and (A) is the safer choice, as it can 
be seen as an expansion in ry ^ 1. Exciting application 
of the extension by means of ri{<f) are found in lower¬ 
dimensional systems [6TH63] . 

The results of this work provide a solid basis for 
studies of competing ordering phenomena in the realm 
of fermionic quantum phase transitions with the ERG. 
Therein, the beta functions for the 0(A^)-|-0(M)-model 
appear as a subset of the larger set of flow equations for 
the whole system. For instance, such a system could be 
given by the Gross-Neveu-Yukawa-model with 0(1)- and 
0(3)-symmetric order parameter fields in 2-1-1 dimen¬ 
sions, describing multicriticality of gapless Dirac fermions 
in graphene [5^. Within LPA' it will be exciting to 
see whether associated scaling relations are still satis¬ 
fied within (A). Furthermore, our finding of the scheme- 
independence of 03 for both the IFF and DFP - if it 
persists in the presence of fermions - allows to unambigu¬ 
ously resolve the corresponding value of 03 and thus the 
stability of multicritical points from the flow. Note that 
fermion-boson-couplings are in general also ((i-dependent, 
e.g. given by a Yukawa coupling which introduces 

new challenges for solving the scaling equations [HU [65] . 
Again, Taylor expansion schemes or spectral methods 
qualify as candidates to address such questions. 


Appendix A: Recursion relation for expansion 
coefficients 


The coefficients Vn for the single-field models can be 
derived from recursive relations. To obtain the expres¬ 
sions we employ Faa di Bruno’s formula 

n 

= Y, (Al) 

fc =0 


Here ... ,Xn-k+i) are the Bell polynomial de¬ 

fined by 



(A2) 


with the sums being over sets {rrij} = {mi ,..., mn-k+i) 
such that mj = k and This combi¬ 

natoric definition, however, is not of practical relevance 
since the Bell polynomials are implemented in many com¬ 
puter algebra system. For instance, in Mathematica 8 
[HH] they can be accessed via 

... ,X„_,.+i) = BellY[n, k, {A,}]. (A3) 

The number of elements in the set {Xj} = {Xi,..., X„} 
can be larger than n — k since only the first n — k elements 
are used for the evaluation of B^^k- We set Bq^q{X) = 1. 

We now derive the Taylor coefficients for the 0(A)- 
model. The recursion is initialized by vq = Nb/d{l + a) 
and V 2 = <J. All coefficients with odd n vanish. The right 
hand side of Eq. ([^ is given by 

/3(0) = - dv{(j)) + a(j)v'{(j)) + bf{v"{(j))) 

+ {N-l)bf{q{^)), (A4) 


where f{x) = (1 -I- x)~^ and q{(j)) = 4>~^v'{4>). With the 
ansatz ( [T^ for v{(l}), the expansion coefficients of v"{4>) 
and q{(p)&X (p = 0 are given by v „+2 and Vn+ 2 l{n + 1), 
respectively. With Eq. (Al I we find that /3„ = ;d("’)(0) is 
given by 


/ 3 „ = {an - d)Vn + fkBnA'^S, ■ ■ ■,Vn+2) 


+ {N-l)bYfkBn,k{Y, 

k=0 


Vn+2 \ 

n -I- 1/ 


(AS) 
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with 


fk = 


{-lfk\ 

(1 -I- V2)^^^ ' 


(A6) 


The coefficients of the scaling solution are now found 
from Pn = Vn = 0. We note that Eq. (AS) is linear in 
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the highest coefficient Vn+ 2 - This is a result of the fact 
that Bn^k{Xi,... ,Xn-k+i) can only depend o n X n for 
k = 1. In this case, however, we find from Eq. (A2) that 
the polynomial equals Consequently, Vn +2 is found 
from 


0 = {an — d)vn + Nhfo — 


bv. 


n+2 


(l + a)2 

n 

+ b'^fk Bn,k{v3,---,Vn+l) 


0- 


TV- 1 
n + 1 


k=2 


(V 3 

u„+iy 

V 2 

n ). 


(A7) 


Appendix B: Flow equations for the two-field model 

We recall the flow equations for the effective poten¬ 
tial and the anomalous dimensions of the 0(A)-|-0(M)- 
model within LPA' for the optimized regulator [iniiii]. 
We write 


p = ^2/2, M = xV2. (Bl) 

The dimensionless effective potential u is defined via 

u{p,fj.) =v{'iIj,x)- (B2) 

The beta function of u is given by 
/3{p, p) = — du + 2aipup + 2 a 2 pup 

_l_ i>i{\ + Up2pUpij) b2{\ -\-Up2pUpp) 

(1 -I- Mp -I- 2pupp){l + Up + 2pupp) - 4:ppu^pp 

bi{N-l) b2 {M-l) 


1 + Up 


1 + Up 


(B3) 


such that u{p,p) = fd{p,p). 

To compute the stability matrix of the two field model 
we apply the expansion 


K K-r 




n—0 m—0 


n!m! 


(p-por(p-por 


(B4) 


where {po,po) is the location of the minimum of the ef¬ 
fective potential. Its scale dependence is given by 


Po = 


Po = 


uo2/3^^’°Hpo:Mo) - Po) 

U 20 U 02 - u^i 

U2oI3^°’^\po, Po) - UiiP^^-^\po, Po) 

U 20 U 02 - u^i 


(B5) 

(B6) 


The flow equations for the coefficients Unm with n,m > 2 
read 


'dnm — ’ '^{pOiPo) T '^n+l,mP0 “f '^n,m+lP0' (B7) 

Again, the implicit dependence on the location of the 
minimum has been taken into account. The set of run¬ 
ning couplings is given by 


{f/nm} — (pO; MO) {Unm}n-|-m>l)- (B8) 

The elements of the stability matrix read 


(A): 


M 


Timin'm' 


{B) : 


M 


Timin'm' 



for the schemes (A) and (B), respectively. The expres¬ 
sions are evaluated for u*. This generalizes Eqs. (27) 
and (28) to the case of two order parameter fields. 

The anomalous dimensions for the two field model read 


16ud pou\^ + po(u2o -k 2/xoA)^ _ 

d (1-1- 2poU2o + 2poUo2 + 4poPoA)2 ’ 

16vd poufi -k po(uo2 -k 2poA)^ 
d {1 + 2 pqU2o + 2pqUo2 + 4:PqPqA)'^ 


with A = U 20 U 02 — uh- At the DFP we have u*pi = 0 
and A* = u*_ 2 ou*,o 2 ) and the expressions generalize to 
the anomalous dimensions of the individual single field 
models. 
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